Modulation of breathers in the three-dimensional nonlinear Gross-Pitaevskii equation. 
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In this paper we present analytical breather solutions of the three-dimensional nonlinear generalized Gross- 
Pitaevskii equation. We use an Ansatz to reduce the three-dimensional equation with space- and time-dependent 
coefficients into an one-dimensional equation with constant coefficients. The key point is to show that both 
the space- and time-dependent coefficients of the nonlinear equation can contribute to modulate the breather 
excitations. We briefly discuss the experimental feasibility of the results in Bose-Einstein condensates. 
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Introduction - Breathers or breathing solutions are nonlinear excitations which concentrate energy in a localized and oscil- 
latory manner. In various physical systems, such as in Josephson junctions lll|2|], charge density wave systems |3], 4-methyl- 
pyridine crystals [4], metallic nanoparticles [5], conjugated polymers |@], micromechanical oscillator arrays 01, antiferromag- 
netic Heisenberg chains and semiconductor quantum wells llQU . the breather excitations play an important role, directly 

affecting the electronic, magnetic, optical, vibrational and transport properties of the systems. 

In the above mentioned studies, one usually considers genuine breathers, i.e., solutions which oscillate in time when the 
nonlinear equation presents constant coefficients (i.e., without modulation). However, in a more realistic scenario the several 
parameters that characterize the physical systems may depend on both space and time, leading to breather solutions that can be 
modulated in space and time. The presence of nonuniform and time-dependent parameters opens interesting perspectives not 
only from the theoretical point of view, for investigation of nonautonomous nonlinear equations, but also from the experimental 
point of view, for the study of the physical properties of the systems. In this context, in a recent work we have considered 
modulation of genuine breather solutions in cigar-shaped Bose-Einstein condensates (BECs) with potential and nonlinearity 
depending on both space and time, in the one-dimensional (ID) case 111 ill . 

The study of BECs of dilute gases of weakly interacting bosons, realized for the first time in 1995 on vapors of rubidium 
|12] and sodium 11311 . constitutes a very interesting scenario to modulate breathers, since they are well described by a three- 
dimensional (3D) Gross-Pitaevskii (GP) equation arising from a mean-field dynamics 11411 . In the BEC context, one finds high 
experimental flexibility to control nonlinearity via Feshbach resonance, and confinement profile via optical lattices and harmonic 
and dipole traps lll5ll . and there we can also investigate the effects of dimensionality reduction on the soliton solution. 

In the case of a strong trapping in two spatial directions, the 3D GP equation reduces to the simpler one-dimensional (ID) 
form, giving rise to the so-called cigar- shaped configuration. The ID GP equation is a nonlinear Schrodinger equation, which 
can also be used to investigate pulse propagation in bulk crystals or optical fibers 11611 . In a former work, however, the search 
for analytical solutions of the ID GP equation with stationary inhomogeneous coefficients has been implemented via similarity 
transformation [17]. More recently, however, the case of space- and time-dependent coefficients were considered for the cubic 
lUih . the cubic-quintic |[l9ll . the quintic lEoll . and also the GP equation in higher dimensions lEHl . 

The similarity transformation was also used to study self- similar optical pulses in competing cubic-quintic nonlinear me- 
dia with distributed coefficients 12211 . nonautonomous matter- wave solitons near the Feshbach resonance 12311 . bright and dark 
solitons in a periodically attractive and repulsive potential with nonlinearities modulated in space and time [24], solitons of two- 
component Bose-Einstein condensates modulated in space and time l25ll . and quantized quasi-two-dimensional Bose-Einstein 
condensates with spatially modulated nonlinearity [ 26] . 

On the other hand, the search for analytical solutions of the 3D GP equation has attracted a lot of attention due to the fact 
that solutions of higher-dimensional GP equation with constant coefficients are usually unstable 12711 . while the nonautonomous 
GP equation can engender conditions which contribute to stabilize the corresponding solutions 12811 . Taking into account this 
possibility, exact solutions to 3D GP equation with varying potential and nonlinearities were proposed in Ref. 12911 . while in Ref. 
13011 the authors studied analytical 3D bright solitons and soliton pairs in BECs with time-space modulation. 

In this paper, our aim is to show that genuine breather solutions can be modulated by the 3D GP equation with space- and 
time-dependent potential, if one includes time-dependent coefficient to describe the cubic nonlinearity. In this way, we extend 
our recent work [ 1 1] to the more realistic 3D case. To this end, we use an Ansatz that changes the 3D GP equation into specific 
ID equation with constant coefficients, which is easier to solve. As a consequence, however, we have to deal with a set of 
coupled equations, to ensure validity of the similarity transformation. Below we present explicit results for three distinct choices 
of potential and nonlinearity. 
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Generalities - We start with the 3D GP equation given by 

Bib 1 

where ip = ^(r, t), r G M 3 , V = (d Xl d yi d z ), and u(r, t) and ^(t) are real functions representing the potential and the cubic 
nonlinearity, respectively. 

Our goal is to find breather solutions which obey the above Eq. (Q]). Using the Ansatz 

1 p = p(t)e i ^mr,t),T(t)}, (2) 
we can transform the above equation into the ID GP equation 

i$ T = -U (( + G\$\ 2 <S>, (3) 

where <J? T = d$/dr, = d 2 <&/d(, 2 , and G is a constant factor. Using © into (Q} leads to (0, for p, r/, £, and r obeying the 
following equations 

n = |VC| 2 , (4a) 
Ct + (Vt?) • (VC) = 0, (4b) 
+ = 0. (4c) 

Here the potential and the nonlinearity assume the form 

v(r,t) = -7 lt -\V7 1 \ 2 , (5) 

and 

m = ^ 

Note that the potential © and the nonlinearity © are functions which in general depends on the real phase 77(1*, t) and the 
amplitude £(r, £). In this way, one can use the Eq. (|4k) to obtain the general form of (, given by 

C(r, *) = ci(t)x + c 2 (t)2/ + c 3 (£)z + c 4 (t), (7) 

where the coefficients Cj (j = 1,2,3,4) are time-dependent functions, obeying the relationship 

T t =c 2 +c 2 +c 2 . (8) 

Now, substituting © into (0J)) leads to 77 in the general form 

7?(r, t) = dix 2 + d 2 y 2 + d 3 z 2 + ^4x7/ + cfexz + 

+ d 7 x + d 8 y + d 9 z + d 10 , (9) 

where the are time-dependent coefficients which obey the equations 
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where dj = dcj/dt. Next, from d?):) and © we obtain 



■ J{d 1 + d 2 + d 3 )dt 



p(t) exp 

Inserting © into © leads to the potential 

v(r,t) = cdix 2 + (jj2y 2 + W3Z 2 + (jj^xy + UJ5XZ 



(id 



+ w 6 y2; + oj 7 x + oj 8 y + ugz + cu 10 , (12) 
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FIG. 1: (Color online) Plots of the 3D breather solution \ip\ 2 for v(r, t) = and g = — 1, in the (x, t) plane, considering y = z = 0. We 
display the cases (a) c± — C2 — C3 = I/a/3, (b) c\ — C2 = C3 = 1, (c) ci = C2 = C3 = 1.5, and (d) a = C2 — C3 = 2. 



where uj = u(t), with 
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Finally, we substitute and (TTTb into © to get 



flr(t) = Gr t exp 



d z ) dt 



(13a) 
(13b) 
(13c) 
(13d) 
(13e) 
(13f) 
(13g) 



(14) 



Here we note that one can get solutions through the Cj coefficients, that is, for specific choices of potential and nonlinearity, 
one can construct the cj functions, or, for specific choices of cj, one can get the corresponding potential and nonlinearity. We 
illustrate the general situation with the examples below. 

Analytical solutions - To demonstrate the power of the method, we start considering a breather solution of the Eq. ©. A 
specific form of the two-soliton breather solution is obtained for G = — 1, which corresponds to the explicit solution Bill 



*(C,t) = 



air/2 



4(cosh(3Q + 3e 4zr cosh(Q)^ 
(cosh(4C) + 4cosh(2C) + 3cos(4r)) ' 



(15) 



The binding potential of the two-soliton is equal to zero, corresponding to the unstable breather solution. Due to the vanishing 
binding potential, there is a splitting of the solution (1131) into two independent solitons at some point in time Il32ll . On the other 
hand, the spatial and temporal modulation of the nonlinearity and the potential allows that we get stable breather solutions, since 
the solution given by Eq. © presents non vanishing binding potential. Indeed, as we have recently shown [11], the modulation 
of the trapping potential gives support to the coexistence of the two-soliton in the breather solution, without splitting. With this 
in mind, in the following we study the modulation of the breather (TT3T) in 3D spatial dimensions. 

Free evolution - Firstly, we consider the free evolution of the solution (TTSb in the BEC, i.e., we take v(r,t) = 0. In this 
way, we can get the free evolution if we set dj = for j = 4, 5, 10 and dj = —dj/2cj for j = 1, 2, 3. So, we will have 
ujj — —Cj/2cj + 3cj 2 /2cj with j = 1, 2, 3. As an example, one gets ujj = setting Cj = l/y/3 (and C4 = 0), corresponding to 

the following choice of the nonlinearity g = — 1, with p = 1. In this case, the Eqs. 0, ([8]), and © change to ( = (x + y + z)/ \/3, 
r = t, and r] = 0, respectively. Thus, we obtain the breather solution of the Eq. © in the form 

4(cosh(3(x + y + z)) + 3e 12it cosh(x + y + z))e 3it/2 



(cosh(4(x + y + z)) + 4 cosh(2(x + y + z)) + 3 cos(12£)) ' 



(16) 
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FIG. 2: (Color online) Plots of \ip\ 2 for v(r, t) — and g — —1, at the spatial origin (0,0,0). We use — cz — l/>/3 to show the time 
behavior for a — 0.3 in dash-dot line (yellow), a = 0.6 in dashed line (green), and ci = 1 in solid line (red). 




FIG. 3: (Color online) Plots of the nonlinearity (a) and potential (b) for the case with a — 1 + 0.5 cos(t) and c^ — c^ — I/a/3- 

In Fig. Q]we depict the breather solution (fT6b in the (x, t) plane, considering y = z = 0, for (a) c\ = C2 = cs = l/\/3, (b) 
ci = C2 = cs = 1, (c) c\ = C2 = cs = 1.5, and (d) c\ = C2 = cs = 2. Similar behavior is obtained in the (?/, t) or (z, t) plane. 
There we clearly see that the oscillatory frequency increases with the increasing of the values of Cj . We further illustrate this 
fact in Fig. [2 displaying \ip\ 2 at the spatial origin, considering c\ = 0.3, 0.6, and 1, respectively, with C2 = C3 = 1/ y/3. 

Harmonic potential - In this second example we consider the case of harmonic potential. Firstly, we look for solutions 
satisfying dj = —dj/2cj for j = 1,2,3 and dj = for j = 4, 5, ...10. As before, adjusting appropriately the values of Cj one 
can obtain the harmonic potential. We will analyze two cases: i) harmonic potential in a single spatial direction, say x, which is 
obtained with c\ = 1 + 0.5 cos(t) and C2 = C3 = 1/ y/3; ii) harmonic potential in the three spatial directions, which is obtained 
with Cj = l + 0.5 cos(t) for j = 1, 2, 3. 

In the case i), let us have r = 1.79t+sin(t) + cos(t) sin(t)/8,C = (1 + 0.5 cos(t))x + y/3y/3 + V3z/3, p = y/l + 0.5 cos(t), 
r] = sin(t)x 2 /4(l + 0.5cos(t)), and v(x,t) = - cos(t)x 2 /4(1 + 0.5cos(t)) - 0.375 smCffx 2 /(l + 0.5cos(») 2 , which is 
physically implemented by a time-dependent harmonic potential plus an optical superlattice 113 311 . The nonlinearity and potential 
are displayed in Fig. [3^ and[3j), respectively. Also, in Fig. [4] we depict the breather solution for the case i, in the plane (a) (x, t), 
(b) (y, t), and (c) (z, t), respectively. Moreover, in Fig. [5^ we plot the breather profile at the origin. 

Next, we consider the case ii). Here we have r = 3.375t + 3sin(t) + 0.375 cos(t) sin(t), ( = (1 + 0.5cos(t))x + (1 + 






(a) 





(b) 





(c) 


16 














1 


»n 
i> 


■ 


m 
t> 




m 
i> 


■1 


in 


1 


+-» in 








8 
















<n 

oi 


• 


»n 

<N* 


• 


IT) 


• 






H 










4 


O 


2-101: 


O 

I 


2-101: 


O 

2 


2-1012 




X 




y 




z 






FIG. 4: (Color online) Plots of the profile of the breather solution in the (a) (x, £), (b) (y, £), and (c) (2, t) plane, in the harmonic potential for 
the case i. Similar behavior appears in the case ii. 
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FIG. 5: (Color online) Profile of the breather solution for the cases (a) i and (b) ii, at the spatial origin (0, 0, 0). 
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FIG. 6: (Color online) Plot of the breather solution in the (x, t) plane, for y = z = 0, in the presence of the linear potential. Similar behavior 
appears in the other (y, t) and (z, t) planes. 

0.5cos(t))?/+(l+0.5cos(t))z,p = (1+0.5 cos(t))(3/2)/2, 77 = sm(t)[x 2 +?/ 2 +z 2 ]/4(l+0.5 cos(t)), g = -3/(1+0.5 cos(t)), 
and v(r, t) = (cos(t) 2 — cos(t) — 3/2) (x 2 + y 2 + z 2 )/ (4 + 4 cos(t) + cos(t) 2 ). This nonlinearity can be obtained experimentally 
in BEC by time modulated Feshbach resonance, but an optical trapping switching from red-detuned to blue-detuned laser-beam, 
and vice-versa, in a periodic fashion will be required for implementation of the time-dependent potential. In Figs. and[5]3 we 
depict the profile of breather solution in the (x, t) plane and at the spatial origin, respectively. 

Linear potential - As a third example, let us consider the case in which the BEC is trapped by a specific linear potential. Here 
we search for solutions of the Eq. (fT2l) with uoj = for j = 1, 2, 6. For simplicity, we also consider uio = 0. To this end, 
we choose dj = for j = 1, 2, 6, dio = 0, c\ — C2 = cs — l/\/3, C4 = sin(t), and dj = —c^/Scj for j = 1, 2, 3 to satisfy 
the Eq. (llOdK In this case, we get r = t, ( = (x + y + z)/y/3 + sin(t), p = 1, 77 = — cos(t)(x + y + z)/y/%, g = —1, and 
v(r,t) = -sin(t)(x + ?/ + z)/V / 3-cos 2 (t). 

In Fig. [6] we depict the breather solution in the (x, t) plane, in the presence of a linear potential. The solution presents similar 
behavior in the other (y, t) and in the (z, t) planes. Here we note that the choice C4 = c±(t) ^ makes the center of mass of the 
solution to move, as expected. 

Ending comments - In this work we have studied the presence of breather solutions in the 3D GP equation, in the case 
of space- and time-dependent potential, with cubic nonlinearity described by time-dependent coefficient. We have obtained 
analytical solutions through an Ansatz which changes the 3D equation into specific ID equation. The results show that the 
breather solution can be nicely modulated in space and time. We have considered three distinct examples of potential and 
nonlinearity: the free evolution of the breather, the case of the presence of harmonic potential and another one, in which the 
system is driven by linear potential. The breather solutions can be controlled through the presence of external apparatus, and 
this may motivate new research in the field since the modulation can generate stable excitations. 
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[1] E. Trias, J. J. Mazo, and T. P. Orlando, Phys. Rev. Lett. 84, 741 (2000). 

[2] P. Binder et al, Phys. Rev. Lett. 84, 745 (2000). 

[3] H. Kleinert and K. Maki, Phys. Rev. B 19, 6238 (1979). 

[4] F. Fillaux, C. J. Carlile, and G. J. Kearley, Phys. Rev. B 58, 11416 (1998). 

[5] H. Portales et al, J. Chem. Phys. 115, 3444 (2001). 

[6] S. Adachi, V. M. Kobryanskii, and T. Kobayashi, Phys. Rev. Lett. 89, 027401 (2002). 

[7] M. Sato et al, Phys. Rev. Lett. 90, 044102 (2003). 

[8] R. Morisaki et al, J. Phys. Soc. Jpn. 76, 063706 (2007). 



6 



[9] E. Orignac et al, Phys. Rev. B 76, 144422 (2007). 
[10] F. Haas etal, Phys. Rev. B 80, 073301 (2009). 

[11] W. B. Cardoso, A. T. Avelar, and D. Bazeia, Phys. Lett. A 374, 2640 (2010). 
[12] M. H. Anderson et al, Science 269, 198 (1995). 
[13] K. B. Davis etal, Phys. Rev. Lett. 75, 3969 (1995). 

[14] L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Clarendon Press, Oxford, 2003). 

[15] C. J. Pethick and H. Smith, Bose-Einstein Condensation in Dilute Gases (Cambridge University Press, Cambridge, 2002). 

[16] Y. S. Kivshar and G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals (Academic Press, New York, 2003); G. P. Agrawal, 

Nonlinear Fiber Optics (Academic Press, New York, 2001). 
[17] J. C. Bronski et al, Phys. Rev. E 63, 036612 (2001); J. C. Bronski et al, Phys. Rev. E 64, 056615 (2001); V. A. Brazhnyi and V. V. 

Konotop, Mod. Phys. Lett. B 18, 627 (2004); J. Belmonte-Beitia etal, Phys. Rev. Lett. 98, 064102 (2007). 
[18] J. Belmonte-Beitia etal, Phys. Rev. Lett. 100, 164102 (2008). 

[19] A. T. Avelar, D. Bazeia, and W. B. Cardoso, Phys. Rev. E 79, 025602(R) (2009); J. Belmonte-Beitia and J. Cuevas, J. Phys. A: Math. 

Theor. 42, 165201 (2009). 
[20] J. Belmonte-Beitia and G. F. Calvo, Phys. Lett. A 373, 448 (2009). 
[21] V. M. Perez-Garcfa, P. Torres, and V. V. Konotop, Physica D 221, 31 (2006). 
[22] J.-F. Zhang et al, Phys. Rev. A 81, 023832 (2010). 

[23] V. N. Serkin, A. Hasegawa, and T. L. Belyaeva, Phys. Rev. A 81, 023610 (2010). 
[24] W. B. Cardoso, A. T. Avelar, and D. Bazeia, Nonlinear Anal. RWA 11, 4269 (2010). 
[25] W. B. Cardoso et al, Phys. Lett. A 374, 2356 (2010). 
[26] D.-S. Wang et al, Phys. Rev. A 81, 025604 (2010). 

[27] C. Sulem and P. Sulem, The Nonlinear Schrbdinger Equation: Self-Focusing and Wave Collapse (Springer, Berlin, 1999). 

[28] H. Saito and M. Ueda, Phys. Rev. Lett. 90, 040403 (2003); F. K. Abdullaev et al, Phys. Rev. A 67, 013605 (2003); M. Matuszewski et 

al, Phys. Rev. Lett. 95, 050403 (2005). 
[29] Z. Y. Yan and V. V. Konotop, Phys. Rev. E 80, 036607 (2009). 
[30] Z. Y. Yan and C. Hang, Phys. Rev. A 80, 063626 (2009). 
[31] J. Satsuma and N. Yajima, Prog. Theor Phys. Suppl. 55, 284 (1974). 
[32] B. A. Malomed, Soliton Management in Periodic Systems (Springer, 2005). 
[33] S. Peil et al, Phys. Rev. A 67, 051603(R) (2003). 



